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ABSTRACT 


Coupled  nonlinear  wave  equations  are  derived  for  the  evolution  of 
the  components  of  the  electric  induction  field  D  in  a  class  of  rigid 
nonlinear  dielectrics  governed  by  the  nonlinear  constitutive  relation 
E  =  X(D)D  ,  where  is  the  electric  field  and  \>0  is  a  scalar-valued 
vector  function.  For  the  special  case  of  an  infinite  one-dimensional  dielec¬ 
tric  rod,  embedded  in  a  perfect  conductor,  it  is  shown  that,  under 
relatively  mild  conditions  on  X  ,  solutions  of  the  corresponding  initial¬ 
boundary  value  problem  for  the  electric  induction  field  can  not  exist 
globally  in  time  ir.  the  L„  sense  if  it  is  assumed  that  the  electric  field 
in  the  rod  is  perpendicular  to  the  axis  of  the  rod  and  varies  as  the 
coordinate  along  that  axis-  Ir  is  also  shown  that,  when  the  initial 
electromagnetic  field  in  the  rod  has  compact  support,  Riemann  Invariant 
arguments  may  be  applied  to  show  that  the  space-time  gradient  of  the  non¬ 
zero  component  of  the  electric  induction  field  must  blow-up  in  finite  time. 
Dome  growth  estimates  for  solutions,  which  are  valid  on  the  maximal  time- 
interval  of  existence  are  also  derived;  these  are  valid  in  the  simple  but 


physically  important  case  where  A(D)  =  A  +  A  IjDJj".  We  also  discuss  re¬ 
lations  with  recent  work  on  the  phenomena  of  self- focusing  and  self-trapping 
for  high  intensity  laser  beams  in  a  dielectric  medium. 
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1.  Evolution  Equations  for  a  Class  of  Nonlinear  Dielectrics 

'theories  of  material  dielectric  behavior  are  based  upon  a  set  of  field 
equations  (Maxwell's  equations)  and  a  set  of  constitutive  relations  which  hold 
among  the  electromagnetic  field  vectors.  In  a  Lorentz  reference  frame  (x^t), 
i  =  1,  2,  3,  where  the  (x" )  represent  rectangular  Cartesian  coordinates,  and  t 
is  the  time  parameter,  the  local  forms  of  Maxwell's  equations  are  given  by 


(1.1) 


3B 

— ~  +  curl  E  =  0 ,  div  P  =  U , 
3D 

curl  H  -  =  0,  div  D  =  0, 


provided  that  the  density  of  free  current,  the  magnetization,  and  the  density  of 
free  charge  all  vanish.  In  (1.1),  Q,  E,  and  H  are,  respectively,  the  magnetic 
flux  density,  electric  field,  and  magnetic  intensity  while  D  =  E  +  P(E), 

>  0  a  physical  constant  and  P  the  polarization  vector,  is  the  electric 

.  3 

induction  field;  the  relations  (1.1)  hold  m  some  bounded  open  domain  ft  c  R  which 
is  filled  with  a  rigid,  nonconducting,  dielectric  substance.  The  precise  nature 
of  the  dielectric  medium  in  ft  is  determined  by  specifying  a  set  of  constitutive 
equations  relating  E,  D,  H,  and  B;  indeed,  without  the  specification  of  additional 
relations  among  the  edectromagnetic  field  vectors,  the  set  of  equations  (1.1) 
represents  an  indeterminate  system. 

fhere  is,  in  existence,  a  wide  variety  of  constitutive  hypotheses  which 
: ave  been  associated  with  theories  of  nonconducting,  rigid,  dielectric  media;  the 
simplest  of  these  is  that  associated  with  the  dielectric  response  of  a  vacuum  in 
which  there  hold  the  classical  constitutive  relations 


2  :  £o 


g. 


where  the  fundamental  physical  constants  p^  satisfy  e0  p^  -  c~‘ ,  o  being 

the  speed  of  light  in  a  vacuum.  In  1873  N&xwell  [1J  proposed  as  a  set  of  con¬ 
stitutive  laws  for  a  linear,  rigid,  statioriary  non-conducting  dielectric  the 
relations 

D  =  e  •  E,  B  =  p  •  H 

where  e ,  p  are  constant  second -order  tensors  which  are  proportional  to  the 
identity  tensor  if  the  material  is  isotropic.  A  set  of  constitutive  relations, 
which  are  still  linear,  but  which  take  into  account  certain  memory  effects  in  the 
dielectric,  were  proposed  by  Maxwell  in  1877  and  subsequently  used  by  Hopkinscn 
[2J  in  connection  with  his  studies  on  the  residual  charge  of  the  Leyden  jar; 
the  Maxwell-Hopkinson  dielectric  is  governed  by  the  set  of  constitutive  relations 
(  j  £  SI) : 

D  (x,t)  =  c  E  (x,t)  +  /t  <J>(t-x)E(x,t)dt 

-  00 

H  =  p-1  B 

where  e  >  0,  p  >  0  and  4>(t),  t  >  0  is  a  continuous  monotonicaliy  decreasing 
function  of  t,  0  <  t  <  «°.  Noting  that  the  Faxwell-Hopkinson  constitutive 
relations  do  not  account  for  the  observed  absorption  and  dispersion  of  electro- 
nugnetic  waves  in  material  non-conductors,  Toupin  and  Pivlin  [31  generalized  the 
constitutive  relations  (1.2)  and  introduced  the  concepts  of  holohedral  isotropic 
and  hemihedrul  dieLectric  response;  while  the  response  incorporated  into  both 
of  these  theories  is  linear,  they  are  more  sophisticated  titan  (1.2)  in  the  sense 
that  nugnetic  memory  effects  and  coupling  of  electric  and  magnetic  effects  is 
built  into  the  constitutive  theory'.  The  qualitative  behavior  o:  the  electric 
induction  field  in  a  rigid  non-conducting  dielectric  exhibiting  holohedral 
isotropic  response  has  been  studied  by  this  author  in  a  series  of  recent  paper's 


(1.2) 


3 


In  this  paper  we  will  be  concerned  with  ini tial-boundary  value  prub Leans 
associated  with  the  evolution  of  the  components  of  the  electric  induction  field 
D  in  a  relatively  simple  class  of  materials  exhibiting  nonlinear  dielectric 
response.  A  rather  general  theory  of  nonlinear  dielectric  behavior  which  allows 
for  both  electric  and  magnetic  memory  effects,  but  still  effects  an  a  priori 
separation  of  electric  and  n agnatic  response,  was  proposed  by  Volterra  [7]  in 
1912  in  the  form  of  the  constitutive  relations 

-f- 

D(x,t)  =  e  •  E  (x,t)  +  V  (E(x,t)),  x  t 

—00 

t 

B(x,t)  =  y  •  H  (x,t)  +  B  (H(x,t)),  x  e  !1 

—00 

The  constitutive  relations  (1.3)  reduce  to  those  considered  in  [2],  L3J  under 
special  assumptions  relative  to  the  functionals  V,  8,  i.e.,  if  B  =  0,  V  is 
linear  and  isotropic,  and  c  -  cl,  y  =  yl,  then  (1.3)  is  easily  seen  to  reduce 
to  (1.2)  ^  the  particular  class  of  nonlinear"  dielectrics  to  be  considered  in  this 
exposition  results  by  specializing  (1.3)  to  the  situation  where  y  =  yl,  y  >  0, 
3=0,  and  electric  field  memory  effects  are  negligible,  i.e., 


(1.3) 


(1.4a) 


D(x,t)  = 

g(£,t)  = 


P(E(x,t) ) ,  x  e  ft 

•V  w7  7 

y  H(x,t) ,  x  e  ft 


V’e  shall  further1  assume 
neighborhood  of  E  =  0, 


that  de  f 


-  i  0, 

JiJ  M 


so  that  in  a  (Euclidean) 


the  relations  (l.ua)  my  lx?  inverted  so  as  to  yield  the 


constitutive  equations 


H(x,t)  =  y']B(2,t)),  x  e  tt 


(1.4b) 
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As  the  vector  function  jF  is  still  completely  arbitrary,  the  constitutive  theory 
defined  by  (1.4b)  is  still  far  too  general  to  provide  a  tractable  system  of 
evolution  equations  for  the  electromagnetic  field  in  ft;  we  will,  therefore, 
confine  our  attention  to  that  special  case  of  (1.4b)  for  which  there  exists  a 
scalar-valued  vector  function  A(g)  such  that  E(£)  =  ACjjpjj,  with  real 
components  £ . .  Thus,  the  final  form  of  the  constitutive  relations  which  define 
the  nonlinear  dielectric  response  to  be  considered  here  is  given  by 


(1.5) 


E(x,t)  =  A(D(x,t))  p(x,t) ,  x  e  ft 

H(x,t)  =  p-1  B(x,t),  x  e  ft  (p  >  0) 


For  now  we  will  simply  assume  that  0  s  A(p  <  v£>  with  A(£)  >  0, 

VP  i  0;  further  assumptions  on  the  constitutive  function  A  will  be  imposed  below. 
Remarks  .It  seems  worthwhile  to  note,  in  passing,  that  electromagnetic  constitutive 
relations  of  the  form  (1.5)  or, to  be  somewhat  more  accurate,  the  inverted  relations 


(1.6) 


D(x,t)  =  €  (E(x,t) )  E(x,t) ,  x  e  ft 
B(x,t)  =  p  H(x,t),  x  e  ft  (p  >  0) 


liave  appeared  in  the  recent  literature;  e.g.,  Rivlin  (SI  considers  (1.6)  and 
indicates  that  in  an  isotropic  material  conforming  to  this  constitutive  hypothesis 
the  dielectric  "constant"  f  must  be  an  even  function  j'  ’he  rragnitude  of  E»  i.e-, 

€=  C(E  •  E) .  Townes,  et.  al.  fit]  considered  the  problem  of  a  high  intensity 
Laser  beam  propagating  through  a  dielectric  medium;  thiy  assume  that  the  high 

intensity  of  the  beam  affects  the  dielectric  "constant"  £  in  such  a  way  that 

'  ) 

the  effective  €  in  the  medium  is  given  by  €=€0+€,  jjEo  where  £^>0  > 

€  > 0  ;  they  then  go  on  to  demonstrate  that  the  presence  of  the  nonlinearity 


may  give  rise  to  an  electromagnetic  beam  which  produces  its  own  wave  guide 
and  thus  propagates  without  spreading  (the  so-called  phenomena  of  self¬ 
trapping  of  the  beam).  Strauss  [lol  and  Whitham  [171  both  consider  a  polar¬ 


ized  wave  with  frequency  go  propagating  in  the  direction  e  (parallel  to 

the  x^-axis  in  our  cartesian  coordinate  system)  in  a  dielectric.  They  assume 

that  the  high  intensity  of  the  electromagnetic  field  in  the  beam,  given  by 

ikxn  .  . 

„  ,  \  ?  -lU't 

E  =  u(x  ,x  )fxv)e  *e  e 

again  produces  changes  in  the  dielectric  constant  6  so  that  6  =  €q  +  ||Ej|‘ 

Using  a  paraxial  approximation,  i.e.,  |u  |  <  <|k  u  |  ,  k=u> 

X,  x  x„  o 

0  j  -5 

these  authors  ([l6],  [ 17 3 )  then  claim  that  u(x^,x0,x„)  satisfies  a  nonlinear 
Schrodinger  equation  of  the  form 


o 


O 

b  u 

o 

bA 


Using  this  last  equation  (and  setting  x?=t)  these  authors  show  that  under 
an  appropriate  set  of  assumptions  the  intensity  |  uj‘  of  the  beam  blows  up  at 


a  finite  value  of  t  and  thus  claims  to  have  a  rigorous  demonstration  of  the 


phenomena  of  self-focusing  of  an  electromagnetic  beam.  We  will  indicate, 


following  the  statement  and  proof  of  our  first  Lemma  below,  why  we  feel  that  the 
reductions  of  the  pertinent  evolution  equations  for  the  electromagnetic 

field  in  the  beam,  to  the  nonlinear  Schrodinger  equation  (given  above)  for  the 
intensity  u(x^,x0,x^)  ,  are  in  error  and  ignore,  in  effect,  the  basic  non¬ 
linear  character  of  the  dielectric  medium  in  which  the  beam  is  propagating. 
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Lemma  1.  Let  fi  c  be  either  a  bounded  or  unbounded  domain  and  assume  that 
>2  is  filled  with  a  rigid,  nonlinear,  nonconducting  dielectric  substance  which 
conforms  to  the  constitutive  hypothesis  (1.5).  Then,  in  ft,  the  components 
D^(x, t) ,  of  the  electric  induction  field,  satisfy  the  coupled  system  of  non¬ 
linear  wave  equations 

(1.7)  m— =  V2(A(D)D.)  -  —  (grad  A(D)  •  D) ,  i  =  1,7, 7. 

3t  i 

Proof.  We  begin  with  the  identity 

A  A  =  grad(div  A)  -  earl  earl  A 

which  is  valid  for  any  sufficiently  smooth  vector  field  on  SI;  applied  to  the 
electric  field  £(*,t)  the  identity  yields 

(1.8)  V2E.  =  tt'—  (div  E)  -  (curl  carl  E) . ;  i  =  1,2,3. 

l  3x.  ~  ~  l 

L 

In  view  of  Maxwell's  equations  (1.1),  and  the  second  constitutive  relation  in 
( L 5 ) ,  we  have 

3B 

curl  curlE  =  -  curl  (~) 

~  d  t 

3H 

-  -  v  curl  (v--) 


-  -U  ~  (curl  H) 


3‘p 


-’j 


31 


2 


so  that  (3.3)  h.as  the  equivalent  form 

3 


(1.9) 


32E>.  ~ 

»—T=  v  h  -  r^7  (div^>  1 

dL  i 


=  1,2,3 


3F  3D  9D,. 

3E^"  '  3xT  Ajk  (2}  Sxf 


,  -  2F.  ,  3D. 

V  E.  =  r-d-  ( ■=■ — - )  =  -i—  (A..  (D)  ~  ), 
l  3x^  dx^.  3  £  ~  3Xy 


A 


<y  (l.'.b), 


7 


J£'l 

where  A.  .(D)  5  ,  and  the  standard  summation  convention  has  been  er.idcyed. 

Thus  (1.9)  becomes 


(1.10) 


32°i  ,  3D- 

(A.  ,  (D)  ~  ) 
1 11  ~  9x, 


,.2  9x,  i£ 

d  K 


a  an 

357  (A]h  3^7  > 


However,  by  virtue  of  our  hypothesis  that  E. (D)  =  A(D)  D. ,  we  easily  find  that 


and  therefoi'e 

(1.11) 


A.  .  (P)  =  A ( D)  6  .  .  +  ~  D. 
n  ~  l"1  3D.  i 

J 


92n. 

t 


9D„ 


~t  -  (  i  A(rj)  6..  t~p.i4) 

3tz  :"k  I?"  d  1 


3\ 


~  (  [  A(D)  6..  +  D.l  ~  ) 
•JX^  ~  -)h  dD^  J  3x. 


where  we  sum  on  each  repeated  index;  expanding  (1.11)  and  using  the  Maxwell 
3D. 

relation  div  D  =  =  0,  we  obtain  the  stated  result  (1.7) ,  i.e. , 

j 

7^ 

^  •  ’x  an.  «x  <x x  r, % 

(1.17)  U - 7  =  J-  ( A ( D )  )  +  y—  ’  (  -r— '  !)•)  -  ~~  (  — •  P.) 

2  9x,  ~  9x,  9x,  3x,  x  9x.  9x.  : 

}.  k  h  i  ] 


9t 


x 

■*» 

9~ 


^ — v —  (A(P)D.)  -  J-  (  ~P.) 

JX,  dx,  ~  1  OX.  OX.  •] 

r.  k  i  ;  J 


Remarks  , , 

We  now  retux’n  to  the  di  scussion  of  '.he  woi  k.  of  Strauss  ll6],  and 
Whithaxn  1 17 j  which  we  began  prior  to  the  statement  of  the  Lemma  above, 
if  we  take  the  constitutive  equations  (l.'>)  in  the  inverted  form  (l.o)  and 
substitute  into  (1.9)  we  obviously  obtain 


M>  — p  (€(E)e7  -  7"~E^  -  $/bx^(div  E) 
bt 

However,  in  a  nonlinear  dielectric  media  it  is  not  generally  true  that  div  E  =  0 
'--d  thus  the  term  b/bx.j(div  E)  can  not  be  discarded  in  the  above  evolution 


equation  for  the  electric  field.  In  particular,  if  D  =  £(E)E  =.  (£  +  )t 

then  div  [  (€  +  CjlEH^E)  ]  =0  and  not  div  E  =  Q  .  In  [  lr;  ]  however  (page  b!*9' 
the  author  has  tacitly  assumed  that  div  F.  =  0  even  though  he  proceeds  to 
employ  a  nonlinear  constitutive  relation  between  P  and  E  (and,  thus,  between 
D  and  E)  while  in  [ 16 ]  the  author  begins  with  the  standard  Electromagnetic 
wave  equation, 


_x 


c 


(€  E)  =v2  E 

at" 


7  -iwt 

assumes  the  form  E  =  u(x  ,x,,,x  )e  -'e  e 
last  equation  to  an  equation  for  u(x  .x^x. 


for  the  wave,  so  as  to  reduce  the 
of  the  form 


f>ik  ^7  -  A  u  +  (k  -  ^|-)u  =  0  , 


a* 


3 


and  then  assumes  that  the  high  intensity  of  the  (laser)  beam  modifies  the 


dielectric  character  of  the  beam  so  that  €  =  €p  +  60  !|Ejp;  this  form  for  £  is 
then  substituted  into  the  last  equation  for  u(x, ,xn,x_ )  so  as  to  give 

1  c  $ 

(modulo  an  approximation)  the  "appropriate"  nonlinear  Schrodinger  equation  for 
u  •  The  problem  with  all  of  this  is  that  at  that  point  at  which  the  beam 
has  modified  the  character  of  the  dielectric  medium  in  which  it  is  propagating, 
so  that 


IN 

o 

it  is  no  longer  true  that  div  E=0  and  the  standard  equation  ,-,(££)  =  7  E 

ikx„  c"  ftt 

is  no  longer  valid,  i.e.,  the  assumed  form  for  E  ,  E  =  ue  ~e~  w  e  must  be 
substituted  into  the  more  general  equation 


M- 


i/5xi  (div  E) 


I 


with  c(E)  |E'  .  Under  these  conditions  the  ochroiinger  equation 

derived  in  [lr],  [  16 j  for  the  intensity  u(x^,xj,x^)  will  clearly  not  result 
and  does  not  seem  to  follow  from  any  reasonable  set  of  approximations ;  a 

ikx. 

i  -IV 

rigorous  demonstration  of  s.:lf-fccusing  for  the  beam  described  oy  l  -  ue  'e 
would  therefore  seem  to  be  an  open  problem. 

We  now  assume  that  cfi  is  sufficiently  smooth  to  admit  of  applications 
of  the  divergence  theorem  and  we  denote  by  v(x)  the  exterior  unit  normal  to 
5Q  at  a  point  ;  we  also  denote  by  t(x)  a  generic  vector  in  the 

tangent  plane  to  ft at  x € c f)  •  -ft6  evolution  equations  (1-7)  are  to  hold 
in  some  cylinder  flx!0,T'i,  T>C  ,  in  R  and  wc-  now  associate  with  this  system 
a  set  of  initial  and  boundary  data.  In  Q  we  require  that 


3D. 


(1.13) 


j.(x,0)  =  f.(x),  -  ( x , 0 )  =  g.(x),  J  x  e  ft 

1  1  ~  dl  ~  1  v 

!  i=l  ,2  ,3 


Standard  results  from  electromagnetic  theory  L 9 ,  §131  also  dictate  that  if  u 

3 

is  a  bounded  domain  in  P  then 

(1.14a,'  [D(x,t)  •  v(x)  i  =  o(x),  (x,t)  e  of)  x  [0,T) 

'"..14b)  [£(x,t)  •  t(x) ]  =  0,  (x,t)  c  31.  x  id, 7) 

In  the  s-»t  c*  relations  (1.14),  I r(x)  !  denotes  fie  jump  of  the  scalar-valued 

f unction  1  across  at  x  c  dft  while  o(x)  denotes  the  density  of  surface 

charge  it  the  point  x  e  31;  these  boundary  conditions  can  be  written  m  an 

-j !  t •  rTat ive  form  as  follows:  If  we  let  b  (x,t)  denote  the  '-i^ctric  induction 

field  at  points  (x,t)  e  P.3/fi  x  [0,T)  then  (l.lha),  (1.14b)  are  clearly 


equivalent  to 
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(1.16a)  D(x,t)  •  v(x)  -  P  (x,t)  •  v(x)  -  a(x),  (x,t  )  c  iw  *! 0,T) 

(1.15b)  X(r(x,r))D(x,t)  •  t(x)  =  lf(x,t)  •  t(x),  (x,t  )  t  3d  x  (0,T) 

where  u  (x,t),  (x,t)  t  R' /’.l  x  0,7),  !■;  the  elc  trie  field  associated  wit; : 

-  (x,t).  In  part icular ,  . :  d  ii  e  p.  ,  ^nd  fi/Sii  is  filled  with  a  perfect 

eon'iuctor  (in  which  h  -  ]i'  -  d)  *.h-.:.  (L.ll,a),  (1.15b)  co<lu<  <•  to 

(,1.16a)  D(x,t)  •  v(x)  -  a(x),  (x,t)  r  >  [0,T) 

a.Kb>  A(D(x,t) )  P  (x,t  )  •  t(x)  =  0,  (x,t>  <  311  x  fO,T). 

In  this  p&pur  we  wic-h  to  consider  tlu-.t  p.ir*<  leular  subcase  of  the  general 
. r.j :  . i .  -boundary  value  problem  (1 ,  ( l .  1  o ) ,  ( 1  .  lui,b)  which  corresponds  to 
tho  assumption  that  the  geometry  of  h  is  an  infinite  one-dimensional  (non¬ 
linear  dielectric  rod);  we  want  to  investigate  whether  a  smooth  electric  field, 
which  is  perpendicular  to  the  axis  of  the  red,  and  depends  only  on  variations  of 
the  coordinate  along  that  axis,  can  exist  globally,  i.e.,  for  t  fc[0,*).  We  assume 
therefore,  that  the  rod  occupies  the  cot.f igurntion  depicted  in  Figure  1,  below. 

The  problem  of  considering  a  finite  rod  glees  rise,  as  a  consequence  of  the 
.n-pronr  i  ate  ••pecia!  i  rat  lops  of  ( t .  1 6 )  .  to  ch>  imposition  of  a  priori  smoothness 
assumptions  on  D(::,t)  at  the  phnar  boundaries  of  the  rr.J.  We  comment  on  this 
situation  at  the  end  of  §3. 
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(1.15a)  D(x,t)  •  Ax)  -  D  (x,t)  •  v(x)  -  o(x) ,  (x,t)  e  0,T) 

(1.15b)  A(P(x,t))D(x,t)  •  t(x)  =  I)  (x,t)  •  t(x),  (x,1)  c  oil  ■<  (u,T) 

where  L  (x,t),  (x,t)  t  R'./W  x  >  0  ,T) ,  ‘lie  electric  field  associate.’.  wit'?; 

D  (x,t).  In  punt  Lculur,  if  d  ■-  i)  c  R' ,  and  ii/Si  is  filled  with  a  j.erfe:  f 

«  :’c 

com juc tor  (in  which  b  -  h  -  d)  t.h  ;>.  (1.15a),  (i.lSb)  redwf<'  to 

(1.16a)  D(x,1  )  •  v(x)  -  a(x),  (x,t)  e  (>S1  x  [0,D 

(1 .161)  A(D(x,t))  P (x,! )  •  £(x)  -  0,  (x,t)  *  fO,T). 

In  this  paper  we  wish  to  consider  that  particular  subcase  of  the  general 
.ni1- 1 i '.-boundary  value  problem  (1.7 ),  (1.1  i),  (:.iGn,b)  whicii  cormyondt  to 
the  assumption  that  the  geometry  o  f  ll  is  an  infinite  one-dimensional  (non¬ 
linear  dielectric  rod);  we  want  to  investigate  whether  a  smooth  electric  field, 
which  is  perpendicular  to  the  axis  of  the  rod,  and  depends  only  on  variations  of 
the  coordinate  along  that  axis,  can  exist  globally,  i.e.,  for  t  ([0,").  We  assume, 
therefore,  that  the  rod  occupies  the  oonf iguration  depicted  in  figure  1,  below, 
ihe  problem  of  considering  a  finite  rod  gives  rise,  as  a  consequence  of  the 
i;  propr ; ate  specializations  of  (1.1b),  to  the  imposition  of  a  priori  smoo ; hue-  , 
assumptions  on  D(x,t)  at  the  planar  boundaries  of  the  i,  '.  We  comment  on  tln- 
situation  at  the  end  of  §3. 
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Specifically,  we  take  for  the  finite  cylinder 

(1.17)  2  =  {(x^x^x^)  I  x.  real,  i  =  1,2,3,  *  ", 

^(x,,x^)  =  C^( const. )} 

with  generators  parallel  to  the  Xj  axis  and  w e-  assume  that  for  soma  small  c  >0 


2  n  {(x  ,x  x,,)  |  -<*•  <  x,  <  <»} 


9  2  2 

C  {(x1,x.,,x0)  I  -»  <  Xj  <  ■»,  x“  +  x  t  L.  ] 


1 


For  ii  vie  then  take  the  (infinite)  circular  cylinder 


(1.18) 


,x  ) 

L  O 


OO  <  <  '*>, 


<5  >  e  >  0} 


and  assume  that  the  annular  region  W/.i  between  the  dieleetrie  rod  and  the 
circular  cylinder  is  filled  witli  a  perfect  conductor;  in  Q  the  dielectric  media 
is  assumed  to  be  governed  by  the  constitutive  hypothesis  (1.5).  We  want  to 
examine  the  possibility  of  there  existing  in  the  rod  a  smooth  electric  field 
which  is  perpendicular  to  the  x^x.,  plane  and  hence,  orthogonal  to  the  axis  of 
the  dielectric;  specifically,  we  are  interested  in  smooth  electric  fields  of  the 
form 


(1.19)  E(x,t)  =  (0,E.,U  ,0.0),  —  ■  x.  . 

Of  course,  in  S2/ 12  we  must  have  E  ~  0.  In  order  to  proceed  with  the  reduction 
of  the  evolution  equations  (1.7),  which  corresponds  to  the  situation  at  hand, 
we  will  need  some  additional  assumptions  relative  to  tin  constitutive  function 
\ ;  specifically,  the  hypotheses  on  \  which  wilt  hold  I hr< 'Ughnw t  the  rest  of 


this  section  are 
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<Xl)  X  €  C1  (R3 ;  [0,-)),  A(C)  >  o,  VC  t  o 

(X2)  | X  ' '  (  c  )  I  <  ®,  Vr,  £  R1 

( X 3 )  0  <  c  X'Cc)  +  X(C)  <  Vc  £  R1 


*  .  1  ft 

where  \\0  ~  X((0,  c,  0)),  r  e  R  ,  By  (XI)  and  the  definition  of  A  it:  is 
immediate  that  X  e  C^CR1;  [0,°°)). 


We  now  proceed  with  the  induction  of  the  nonlinear  evolution  equations 
(1.7).  In  view  of  (1,5^),  (1.19),  in  ft 

(0,C20))  =  A(D)(Dj,D2  ,l>3) 

from  which  it  follows  that,  in  ft,  D1  =  D3  =  0  and  E^x^t)  =  X(D)D2(xx  ,x?  ,x3  ,t) . 


3D 

However,  div  D  =  ~ n  ^  ^ 

~  3:<2  -  0  so  th.it ,  for  each  t  £  0, 

x-j  , x,, .  As  E,;  depends  only  on  x^ 


can  depend,  at  most,  on 


?  A  /  \ 

__  =  ___  (A(D)D2(x],x3,t)) 

~  g~~  (  X  ( 0 ,  ,  C )  Dj  ( x.  ,  x  ^ ,  t )  ) 

=  5x^'  (^(D2(xi’xvt:')D2(xl’x3,t)) 


3D?  ^ 

f  A  A 

k X ’  (D„)D  +  X(l'  )) 

-  0. 

5x^ 

/  2  2  ' 

3D2 

By  hypothesis  (X3)  it 

then  follows  that 

=  0  and ,  thus ,  i  n 

(1.70)  D(x,t)  =  (0,  D^(x^ ,t) ,  0) 


In  view  of  (1.20),  not  only  is  div  D  =  0  automatically  satisfied  in  ft,  but, 

At  this  point  hypothesis  (A3)  could  be  weakened  to  the  assumption  that 
*  1 

(l,\  (O)'  t  0  a.e.  on  R  and  (1.20)  would  still  obtain. 


lh 


as  is  easily  ver,lfied5  so  are  trhe  nonlinear  evolution  equations  \±./)  iv?r 

i  =  1,3,  i.e. , 


(grad  A(D)  •  D)  -  v—  (  ^  MD^x^t))  •  D^x-pt)) 

1  i  2 

=  0,  i  =  1,2,3 

while  (D.)tt  =  V2(A(D)D. )  =  U  for  i  =  1,3.  For  i  =  2  we  then  obtain,  for 
-«■  <  <  <»,  and  0  t  t  <  T, 

(1.21)  y  ^  (Xxt)  =  V2[X(D2(x1,t))  D2(x1,t)] 

32  * 

=  — 0  i'A(D,(jc,  ,t)>D  (x  ,t)  3. 

°X! 

In  view  of  our  assumption  that  the  rod  is  infinite  in  extent,  the  boundary 

conditions  (1.16)  do  not  come  into  plav  here.  In  fact  to  simplify  the  analysis 

3D2 

we  will  now  assume  that  the  initial  data  D^x^.O)  'ind  ^  ‘(x^,11),  °°  ■ 

have  compact  support  on  1?  .  Then,  in  view  of  the  fact  that  hypothesis  (13) 
implies  the  strict  hyperbol i ci ty  of  (1.21),  D^tx^.t)  will  also  have  compact 
support  on  R1,  say  supp  l>2  t  for  some  6,  0  <  C  <  «,  for  as  long  as 

a  smooth  solution  of  the  initial-value  problem  for  (1.21)  exists. 

We  now  set  x,  £  x,  D.,  £  i.  The::,  for  th*-  physical  si  tout  ion  lesc’riiv  1 
above,  the  initial -boundary  value  problem  associated  w: *  li  '  1  e  coupled  system  of 
nonlinear  evolution  equations  (1.7)  reduces  to  the  following  nonlinear,  one- 
i unensiona! ,  initial  value  problem  on  th,-  x  axis:  find  u  -  u(x,t), 

-uo  «■  x  <-  :)  5  t  '  T,  such  that 


,2  .2 

d  n  _  d 

at2  ax' 


-  - X-  (u\(u),  (x,t)  £  (-■ «,®)  x  i'0,D 


(1.22) {  u(x,0)  =  Uq(x),  u+(x,0)  -vn(x), 


O' 


—  on 


X  < 


5 


★ 

where  u  >  0,  A  satisfies  tiie  hypotheses  (XI)  -  (A3),  and  u^(x),  v  (x) 
have  compact  support  on  R* .  Any  smooth  solution  of  (1.22)  will  also  have 
compact  support  on  R^  and,  as  indicated  above,  we  will  take  the  support  of 
u  to  be  in  the  interval  (-d,^)  where  0  <  6  < 
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?.  Global  Honexistence  of  Electric  Induction  Fields 

In  thi:.  section  we  will  demonstrate  tliat  under1  the  »dd'  tiuml  iiyix.tti.-nji 

*  (2) 

on  the  constitutive  function  a(£)  =  A((f',t,0)), 

( A  4 )  For  all  £  e  R ^  and  some  u  >  2 
ft  /p  pA  (p)  dp  2  C^’A(t') , 

smooth  global  solutions  of  (1.22),  i.e.,  solutions  of  (1.2?)  on  (-«*,<»)  x  [0,T>, 
for  all  T  >  0,  will  not,  in  general,  exist;  in  fact,we  will  show  tliat  under 
relatively  mild  assumpti  s  on  the  initial  data, the  L?(-»,a)  norm  of  u(x,t) 
must  be  bounded  from  below  by  a  real-valued  nonnegative  function  of  t  which 
becomes  infinite  as  t  -  t  <  <*>.  Some  growth  estimates  for  solutions  of  the 
initial-boundary  problem  (1.24),  which  are  valid  on  the  maximal  time- 

interval  of  existence,  will  also  be  derived.  In  §3,  we  show  that  under  stronger 

■i. 

assumptions  dti  A^),  than  that  represented  by  (A4),  it  is  possible  to 

demonstrate  via  a  Riemann  Invariant  argument  that  smooth  solutions  of  (1.22) 

cannot  exist  globally  due  to  finite-time  breakdown  of  the  space  time  gradient 

fu  (x,t),  u  (x , t) ) . 
y  t 

before  proceeding  with  the  art: lysis,  let  us  note  that,  if  we  set 
P(?)  =  CA.'C),  4  «  R1,  and  Ur,)  =  jj*  sj>(p)dp  then  tV  <  )  -  >.’((;)  and  hypothesis 

(X4)  is  equivalent  to 

(Vr)  for  all  (  e  and  some  a  >  2  al(b)  a  4  E'vO. 

i:.e  proof  of  the  global  nonexistence  results  referred  to  at  ore  now  proceeds 
via  a  series  of  lemmas,  the  first  of  which  is  just  an  energy  conservation  theorem 
for  the  solutions  of  (1.24).  Thus,  let  S,  0  <  7i  be  any  constant  such  that 

7  '•  ~~  * 

Will  ic  this  hypothesis  Is  satisfied  by  >  ( f )  --  const.  none  ol  our  results  apply 
to  the  linear  wave  equation,  i.e  ,  see  the  footnote  following  Theorem  III. 
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6  >  6  where,  by  assumption,  supp  u(x,t)  -  (-5,6)1  tnen  we  have 
Lerma  2.  If  we  def  ine  the  total  energy  E(t)  of  the  system  (  <  .22)  by 

(2.1)  E(t)  =  (  £xwut(y,t)dy)2dx  +  J*  (/g(x,t)  pX(p)dp)  dx  , 

then  for  as  long  as  smooth  solutions  of  (1.24)  exist, 

i  '»  r  un(x)  * 

(2.2)  Eft )  =  J  f  '  ( f4<ai  v.  (y)dy)"dx  +  /  "  (fQ  pA(p)dp)  dx 

Proof.  In  view  of  the  definitions  of  !/)(?),  £(£), 

(2.3)  E(t)  =  y  (/x  u  (y,t)dy)2dx  +  f6  Z  (u(x,t))  dx 

*  -a.  _uo  X 

Therefore, 

.  j  x 

(2.4)  E(t)  =  uf  (/^cout(y,t)dy)  (  /..u^y/Ody)  dx 

—  OO 

+  /‘5  Z'(u(x,t))  ut(x,t)  dx 

—  OO 

r\ 

=  / 6  (  /x  u  (y,t)dy)  (/xroA^(u(y,t))dy)  dx 

"  V 

+  /'  E’(u(x,t))  u  (x,t)  dx 

—  uO  ^ 

=  /J0  (/^a,ut(y,t)dy)  *(u(x,t)),x  dx 
+■  / 5  E’(u(x,t))  u..(x,t)  dx 

_  on  ‘ 

where  we  have  used  (1.22)  and  the  compact  support  of  u(x,t)  on  R  ,  i.e. 

2 

/x  -Xr  ip(u(y,t) )dy  =  if>(u(y,t)),y  i\. 

'  9y 

=  <|)(u(x,t)  )  ,X 

-  lira  (if>(u(v,t  )),y  I  ) 

p-v_  ^  P 

P-'J 


tp(u(x,t))  -X 


1R 

as  if»(0)  =  0  by  virtue  of  (XI)  and  the  definition  of  ip  .  Therefor1*: 

(2.S)  E(t)  -  /5  Lij/(u(x,t))  fx  m  u  (y,t)]dx 

-  /**  4Ku(x,t))  u,_(x,t)  dx  •»  J*5  I'(u(x,t))  ut(x,t)  dx  =  0 

—  oo  _  oo  *- 

as  v(;)  =  a'(c).  VC  f  R ^ ,  by  definition,  y(0)  =  0,  and  supp  u  1  (-6,5), 

6  <  6.  Equation  (2.2)  then  follows  by  integration  over  L0,t),  the  definition 
of  E(t) ,  and  the  initial  conditions.  Q.E.D. 

Our  next  iemna  is  concerned  with  establishing  a  certain  lit ferent ial 
inequality  for  a  real-valued  nonnegative  functional  defined  on  solutions  u(x,r) 
of  the  initial-boundary  value  problem  (1>22);  namely,  we  have 

Lenina  3.  Let  u(x,t),  (x,t)  e  (-“,«)  x  [0,T)  be  a  smootli  solution  of  (1.24) 
and  define 


(2.6)  F(t)  =  yjT^  (  /^a)u(y,t)dy)2  dx  +  t$(t+to)2 

where  ft,  t  *  0.  If  A(r,)  satisfies  (XI)  -  (X4),  then  for  0  £  t  <  T 
7.1)  F'"  -  (Y+i>r2  -  -2(  ?Y+1 )  F  (6  *  2E(0)) 


--  >  0  ( v/itf 


n!  constant  -jhi.cn  i risen  :ir  \l-  •  constitutive 


assumption  (X4))  and  E(0),  the  initial  energy,  is  given  by  the  light -hand  side 


of  (2.2). 


I 'roof . 


by  ■  irect  diff<  r^ntiation  we  have 
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(?.8)  r’(t)  =  2y  /_i  <Jt  „  u(y  ,t  )dy)  ( f*  ^  u  t  ( y  ,  t  )  Jy  >  .lx  + 

and 

(2.9)  F"  (t )  =  2y  /_t,  (/^rout( yt)dy)/  dx 

♦  Lt,  (f- oo  u(y,t)dy)  (/*  m  utt(y,t  )dy)  dx  +  26. 

Again,  in  view  of  (1.24  ),  the  definition  of  iji(  C  )  ,  £  t  R1  ,  and  the 

compact  support  of  u,  we  have 

(2.10)  F"  ( t )  =  ?u  it  (/^oout(y  ,t)dy)2  dx 

+  2fS  (/x  u(y  ,  t  )dy )  i|>  (u(  x  ,  t )  )  ,  dx  +  26 

s  (/too  ut(y»t)dy)"  dx 

+  2 fS  ~  {(/x  c,  u  (y ,  t  )dy )  i|Ku(x,t)  )}dx 

-  2/6  u(x,t)  i|>(u(x,t))dx  +  26 
=  2p  If  (J*w  ut(y,t)dy)2  dx 

-  2/^  u(x,t)  E  (u(x,t))dx  +  26 

Bv  addinv  and  subtracting  2u f(  E  (u(x,t))  dx  on  the  right-hand  side  of  the 

*  —a. 

last  line,  in  (2.10)  we  obtain 

(,2.11)  F"(t)  -  2 u  J®  (fX  u  v'v,t)dy)-:dx  -  2a  JS  E  Cu(x,t))  dx 

'  -a.  •  -  00  t  ' 

+  ?/Ma E(u(x,t))  -  u(x,t)  E  '  (u<x ,t))  dx  «■  26 

•  *  a)' 

2  2y  it,  (/t.  ut(y,t)dy)2  dx  -  2«J_t,£  (u(x,t)) 
where  we  have  vised  the  hypothesis  (X4)  in  the  form  given  by  (A4).  However,  in 


r 
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v;‘M  of  the  definitions  of  F(t) ,  i.e.  and  ECO ,  ;  t  P  ,  -}«•  , 

in  (2.11)  may  be  replaced  by 

(2.12)  r"(t)  a  2y  J6  (j*M  u+(y,t)dy)?  dx 

-  2a  [ E( t )  -  iJ-  J5  (J*  u  (y,t)dy)?dxJ  +  20 
=  ( 2+a)y  \h  (/x  ,  n.(y,t)dv)2  dx 

-  2a  EO.J;  +  2(1 

where  we  have  used  the  energy  conservation  result  of  Lemrru  2.  Kindly,  we  rewrite 


the  last  i 

.nequality  in 

(2.12) 

in  the 

fom 

(?.13) 

F"(t) 

>  (2+a)  Cu  / 

6  (/x 

CO  J  -  «> 

U  (V2 

L 

t)dVr  dx  +  3) 

-  a  [0 

x  2E(0) ] 

Combining 

(2.3), 

(2.13) 

and  ( : 

2 . 6 )  we 

now  c! 

it  ;in 

(2.14) 

FF"  - 

(a+2) p ' 2 
'  4  >* 

2 

(2+a) 

[-£ 

(  (X 
*  _  <X> 

u(y  , 

i  )  ty ) '  c’.x  +  0(1 

X 

bit 

« v 

(j- 
*  _  >« 

,*  r  v 

+.))'•) 

'■-*  *  -1 

- 

aide 

+  1 1  ( 

') )  > 

- 

( 2  +a ) 

r  /  f  x 

["d'C. 

u  ( y , 

+  t .  i v  )  (  fX  :  1 

"  >  *  '  -  x  r 

+ 

B  ( t  + 

y 

t  )21 

J 


:  )  d  v )  d  x 


1 
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:2+a)  |(  [y/^  (/^  ia(y,t)dyr  dx 


+  3  ( t  f  i . )  j 


*  ^/l(r.  ut(y,t)Jy)‘  dx  +  3]  ) 

~  (/!. oc  c/f.  U  (y,t)dy)  dx 


+  *(t  +  V)7 


a  r(3  +  2E(0)). 


By  virtue,  of  the  Cauchy -Schwarz  inequality  the  {  }  expression  in  the  Iasi 

inequality  in  (2.14)  is  nonnegalive  for  all  t,  0  <  t  <  T,  and, therefore, 

(2.15)  FFT'  -  (-*— )  r'?  >  -  aF(g  +  2E(0)),  0  <  t  <  T 

The  required  result,  i.e.,  (2.7)  now  follows  directly  from  (2.15)  if  we  set 


Y  =  (ct-2)/4. 


Q.E.D. 


Global  nonexistence  of  solutions  to  the  initial -boundary  value  problem 
(1.22)  can  now  easily  be  shown,  to  be  a  consequence  of  the  differential  inequality 
(2.7)  under  various  assumptions  on  the  initial  energy  E(0)  and  the  initial 
lata  ur(x) ,  v'q(x).  To  simplify  the  discussion  v-  introduce  the  notation 

(,2.15c.''  I(uf.)  '  V-f'  (fx  u  (v)dv)‘  Jx 
n  '  _oo  * oo  ( 

(2.16b)  J(un,v_)  -  2 (fx  u,(y)dy)  (f  "  v.(y)d')  -l* 

(tar  first  series  of  results  (Theorems  I  and  ll)  concern  situations  in  which 
e'O)  <0  ;  while  these  results  may  be  of  some  interest  for  their  own  sake  they 
are  not  relevant  to  the  example  in  whirh> paralleling  the  assumptions  in 

Townes,  et.  al.  ,  [it],  Strauss  I’lbl,  ana  •••■hi  them  •  If  ■  we  have 


*U)  =  +  x  ;  KQ  r  0  .  A. ,  c  o 
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(we  must  have  f(0)  ->  0  tor  any  such  A);  results  concerning  the  c'evelopnieni 
ot  singularities  in  the  gradient  (u  ,u  )  for  this  situ.iti in  vouch  ,  i  => 
not  genuinely  nonlinear,  are  presented  in  S3  as  an  app’ ration  <.f  recent 

work  of  Klaineriaan  and  Maida  121]. 

"n-aorvju  I.  L*i*.  u(x,i)  N.  -  ,v.  ion  c;  .  r  .  •  *  .•  '■ 

function  X(0  =  X(  (0,^,0)  sat  isi  ies  (Xi)  -  (X1*).  I:  7«  a.,,v  >  an: 

t  u  (x)  r 

(2.17)  }_jf°  *(p)dp)  dxs-  j  v3(y bb.)'  :.x 

then  Li  k-  =  >:(u;f)  >  0  arid  t  <  sue:  '.hut 

3  uo 

0.28)  ||  u(t)||  2  >  kG(m,  s  <  4  <  t„_ 

L  (-*  ,  O 


where  LO.t  )  denotes  the  maximal  interval  of  existence  <>t  n(x,t  )  and 
■  .tax 

iir.  G(t)  =  i  ». 
r-»-t 

00 


Sertark :  As  a  consequence  of  >. 2.19),  and  the  fact  that  G(t)  tends  ’o  +  ® 
as  t  •*  t  ,  it  follows  t'-at  t  <  t  ,  i.e.,  that  the  maximal  interval  of 

existence  of  u(x,t)  is  finite;  it  nv.iy  happen,  however,  th.it  t  <  tw 
and  thus,  without  stronger  ascumpt  ions  on  A\C> »  we  cannot  conclude  t  licit 
u(x,t'» ,  o:  i';>  derivaf : ven  ,;d  /w-up  ->■:<  r_  .  topr  L.a'a.  :u:t;  .r  t  -*•  ^  'See 

Sail  i  !Dj,  [11  ;  for  <:  relevant,  discus:  ion  of  the  red. a 4  .snip  between  global 
nonexist  aict:  and  f in.it e-time  Lbov-.n.  theoremr.  for  i-oiut-o.c  of  nonlinear 
evolution  equations), 

.-hoof  (Theorem  1).  In  view  of  (2.17),  f(0)  s  0.  Time,  if  we  set  3  =  0  in 
0.7)  this  inequality  reduces  to 

(2.13)  Fn<  t )  f  "  (t)  -  (y+1)  ir/o)  >  0,  0  <  t  <  t  , 

J  O  U  iJkiX 

F  (-)  -  vF  (f\  u(y,t)dyO  dx.  But  (2.1G  is  easily  seen  to  he  equivalent  to 
(- 

'  0  ( t  >  s  C ,  M  t  s  t 


(2.20) 


rrux 


;V*o  successive  j  r it* vi u  t  ions  <a  ( 2  , ,}.  0  )  yie 1 «. ! 


c?.2i)  r  ~Y(t)  <  -y  r  "Y-i  (0)  r'(o)t  +  r~Y  (0),  o  <  t  < 


or,  is  y  >  o,  r  (t)  >  n 


(2.22)  F„(t)  >  - 


,  /-0  ;|J,\ 


T7T\y  Y  =  r.(t),  0  <  t  <  r 
( t'  )\  ’  max 


Clearly,  lim  G(t)  =  +  «>  where 
t-rt 


i/ro(0)\ 

=  Y \fjW'J 


i  . 

Y  J(..0,vcl 


Also,  as  supp  u  c  (-6,6)  and  6  >  6 


(2,23)  t(fL  u(y,t)dy )'  dx  =  f  (/jKy.tJdy)7  dx 


J"_  (x+6)  (j  _  u“(v,t)dy)  dx 
-6  -o 


C  O  2  f  r  0  '"*1 

(j°  (x+fl)"dx):  (/'  ( fX_  u  '  (y , t )dv) "  dx)  5 

-3  -I  ‘ -  5 

-  ■?/;'  -62  2  s 

(26)  (  .  !  j  u  ( y ,  t )  i  v )  d  x ) 

’-6 


(26)  '  2  •  v'26  j  i*fv,t)dy  -  46"  (y,  t)  dy 


and  therefore. 


(2.24) 


!iu  ( t )  ;■ 2  , _ Y  ~-T  In  (j*  luy.thlv)'  dx. 

l/(-6,6)  43^ 


The  growth  estimate  (2.18),  valid  for  0  *'  t  r  t  ,  with  k  -  ,,  now 

max  4 1 1 3 “ 

fellows  directly  from  (2.22),  (2.24),  the  definition  of  rt>(t),  the  fact  that 


24 


u(x,t)  =0,  <  x  <  -0,  and  the  observation  that  (2.17),  1.2.21,'  and,  r.hus, 

(2.18)  also  hold  as  we  let  6  -*•  6.  Therefore,  if  supp  u  (-6,4)  then  for 

0  t  <  t  ;|u(t)'^  2  k(p;6)G(t)  with  <l(t)  given  bv  (2.72).  Q.F.!'. 

i.  (-“,4) 

There  are  several  otlier  rituationa  in  which  tly  r.  m.  ha rone  I  u:  b  n.  a:, 
that  expressed  by  Theorem  I,  follows;  we  will  examine  two  such  set s  of  circum¬ 
stances  below  which  corresponi  to  situations  in  which  we  have,  respectively , 
J(Uq,Vq)  =  f'  and  J(u  ,vf()  <  0,  with  E(0)  <  0  in  both  cases.  Suppose,  first  u 
all,  that  E(0)  <  0  with  v^(x)  s  0,  <  x <  5;  in  this  case  we  nay  choose  p 

such  that  ?E( 0)  +  8  -  0  und  therefore  (2.7)  reduces  to  (2.1°)  with  F  (O 

0  '  0 

replaced  by 

F(t;  t>0,tQ)  =  u/L(/'-  u(y ,t)uy)~  dx  +  Sr(t+t„r 
0  -  <w  ' 


Therefore,  F(t:  80,tg)  satisfies,  for  0  <  t  £  tmov  , 


(2.25)  F(t ;  60,tc>  2 


r  (0;  8Q,t0) 


,  /  r'‘"'WY 


nvax 

1 

Y  5  H(t)  , 


so  that  lim  H(t)  -  +  »  where 
t-n  (tn) 

uc  Q 

F(  n  •  8  t  i 

(2.7b)  t  (t  )  =  ± _ _ _ u 

-ocPV  >  ft  FfTTgfT 

T  -'’pc»'0' 


l/T(u0)  +  Vo 

Y  V  ’  23oto 


We  note  that,  in  view  of  cur  hypothesis , 


—  u,(x) 

(2.?7)  (Jq  ty(p)dp)  dx  >0,  5-6 


is  net  diff  icult  to  sir'W  that  the  minim..:,  value  of  t  ( ;  ,0  is  achieved  at 

U 


“n  wn 


* 


and  that 


(2.28) 


>(V  ~  'J 


/I(uJ 


S„ 


Choosing  in  (2.25)  we  have,  therefore. 


(2.29) 


u/q  (f*co  u(y,t)dy)  2  dx  +  eo(t+tQ)2 


> 


(T(U  )  4-  B  t 
0  10 


1  - 


t  (tn) 

M  0 


--  f 


1 

Y 


for  Ostst 

max 


I(u„.) 


In  view  of  (2.2*0,  (2.28)  we  then  have  the  growth  estimate  (let 


6—5) 


,j  and  i(uQ)  by  (2.16a),  with  6  in  place  of  5.  The  estimate  (2.30)  establishes 

.'.'ohai  nonexistence  of  solutions  to  the  initial-boundary  vain  -  pt>“ien>,  under  the 

hypotheses  (XI)  -  (X4),  for  the  case  where  the  initial  Hita  satisfy 

x  u  (x) 

vq(xj  :  0,  -»<  <  x  <  6  and  ( JQ  pf(t')d.,)  dx  <■  0. 


* 


jp.d  t  <  U )  <.  (j 


Having  examined  the  cases  where  E(L>  £  0  win.  ifn  v( (j; 

V 


with  J(u  ,v» )  -  0,  v-dx!  =  0,  -<*.  <  x  <3  .  we  now  want  to  look  at  the  s 
o’  0  ’  n  ’ 


I  fU.itiO 


E(0)  <  0,  i.e. , 


(2.31) 


/-a.^Oo(X)^(p)dp)  dx  <  “  j  flv'-f.oo  v0(y)dy)‘  dx 


u  rd  ,  tx 


and  J(uQ,vc)  <  0.  In  this  case  we  my  again  choose  8  -  E(j  such  that 


2 E ( 0 )  +  Bn  =  0  ,  so  that  F ( t ;  B^t^)  satisfies  (2.25),  with  8n  -*•  ?L  ,  for 


0  u 


C  <  t  <  t  .  We  note  that  we  now  have 
max 


(2 . ?2) 


t  (t  )  = 

CO  U 


i ( ao )  +  hlo 


2BQt0  -  |J(ug,v0)| 


where 

(2.33) 


0  u.(x) 

Bn  -  -  uf_  (f^oo  vn(y)dyr  dx  "  ?/_£/o  p  ^(P^dD^  dx  >  0 


and  thus  we  must  choose  t  £  t^  where 


(  2 . 34  ) 


'»>it''<Uo,vd! 


0 


♦  i:.  )  relatively  simple  matter  to  show  that  tn(t^)  achieves  a  minimum  at 


tp  *  tQ  -  •  —  (1J(UC  ,VQ)|  4  /-(Ur,  V(,3  t  (U0)  ) 


28 


w»:  d<'note 


t:  i>  )  -  t  then  we  have  the  estuir 

CO  fl  CL 


(2.  3tF 


P/H/L  u<ya)dy)‘di  *  s  ,:*  *  fu)‘ 

♦  S-cVll  J<»„> 


r<-’<Wo  >'it .  i(v  , 

“TTT-F —  - 


n  whc-iv 


1 
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Dr  0  itit  and  the  companion  estimate 
1T&X  * 


(2-377  i  ||u(t)||2  +  B0(t+t  D 

K  !,(-«,<$)  U  U 


I(U0) 


for  Octet  <  t  and  global  nonexistence  of  solutions  to  the  initial- 
max  00  6 

boundary  value  problem  (1.22)  follows  as  in  the  previous  cases.  We  nay  summarise 
the  two  results  corresponding  to  the  situation  where  E(0)  <  0  as 


'Theorem  II.  Let  u(x,t)  be  a  solution  of  (1.22)  and  assume  tluvt:  the  constitutive 
function  A(p)  =  A((0,p,0))  satisfies  (Al)  -  (A4).  Then 


then 

where 


(i)  If  Vp(x)  =  0,  -"<x<£ 

u(x,t)  satisfies,  for  Oct 
is  given  by  (2.27)  with 


a  un(x)  * 

and  f_J L  pX(p)dp)  dx  <  0  , 


<  t  c  /  — -  ,  the  growth  estimate 
itbx  y  ’ 

A  in  place  of  6. 


(2.30) 


(ii)  If  the  initial  data  (u^(x),  v^(x))  satisfy  (2.31)  and 

ll(fM  uQ(y)dy)  (/*«,  vQ(y)dy)  dx  <  0  , 

•■nen  u(x,t)  satisfies,  for  0  <  t  c  t  c  t^  ,  the  growth  estimate  (2././), 
where  is  given  by  (2.23),  t(  by  (2.35),  (2.  :6b),  -.me  \m  =  t^tg)  vJiere 

t  (t^)  is  given  by  (2.32).  In  both  cases  (i)  and  (ii)  .ibove  the  respective 
estimates  (2.30),  (2.37)  imply  that  solutions  of  (1.24)  cannot  exist  globally, 
;  .e. ,  for  t  e  L0,«>) . 


We  now  want  to  consider  situations  in  which 
6(0)  >0  (i.e.,  \(c)  Aq  +  \^2  ,  Aq>0  ,  V,  > 0  ) . 
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Theorem  III.  Let  X(£)  satisfy  (Xl)  -  (X3^  and  define  J1uqJ 

in  Theorem  I.  If  g  (0)  >0  with  fl(u  ,  v  )  >0  and 

-  9  o  o  — - 

/(u  ,V  )  ,,, 

>8e  fo)(3; 

J?(u  ) 


Jlu  ,  V  ' 
<7  0  0 


then  no  smooth  solution  of  (1.2M  can  exist  for  all  te  [0,»). 

Proof-  Assume  that  a  smooth  solution  does  exist  on  [0,®)  •  Then  (2.7) 
holds  Vt  and  we  rewrite  it  as  (set  0  = 0) 


(2.?8)  FF"  -  (y  +  1)F,<?  >-  ?v‘(2y  +  1)F 

o 

where  v‘=2g(0)>0 

As  F’(0)  =J?(u0»vo)  >0  we  have 

(f"y)'(o)  =-y  f_v/y+1)(o)f'(o)  <0 


By  continuity  (F_Y)'(t)<0  for  t.  sufficiently  small.  If(FY)’(t)^0  for 

* 

as  long  as  smooth  solutions  exist  then  3t  = t  such  that 
F~y'  (t)  <0  ,  t<t*  but  (F“y)  '(t')  =0  . 


We  will  show  that  this  can  not  happen.  Since  F(t'  >0  ,  te  fO,t  ]  may  rewrite 


(2.39}  if  Y)”  <'-Yv“  (2Y  +  1)F  '  <'+‘L'1  ;  te  [0,t  ! 

On  [ 0, t* )  .  (F~y) ' (t) < 0  .  Multiply  (0.39)  thru  by  * ( F_Y) * ( t )  ,  te[0,t*] 

to  oil  aln 

’(F_Y)”(F"Y)  '  >hy  v  '-(?y  +  l)F‘^Y+l)  (F~Y)  ' 


This  condition  is  easily  seen  to  require  tlar  _  (j  «X(t)dp)dx  be 

sufficiently  negative;  thus;  neither  this  result,  nr  -inv  ot  our  other  results 
of  a  similar  nature,  apply  to  the  linen’  wave  equation  obtained  by  taking 
>*  (i  )  =  ron.d  - 


* 


29 


T 


JO 


or 


(2.1*7)  -  Y  F' 


'^y+1^f'<-  (v  +  -'.y F-(->+1))1/^ 

—  O 


But  this  clearly  implies  that 


(2.46)  F'  (t)  >  (4Y<?F(t)  ■"Y+1^(t))1//; 


for  as  long  as  smooth  solutions  exist.  Hence 
F(t) 


(2.49) 


dG 


JF(0)(^v2G  +U  y'^G2^1'1  )1/~ 
o 


>  t 


which  inqplies  a  finite  time  of  existence  for  any  smooth  solution  since  the 
integral  on  the  left  hand  side  of  (2.49)  is  convergent  VF 

Gur  last  result  is  a  growth  estimate  for  smooth  solutions  of  (1.24)  which 
is  valid  on  ] 0»tmav  )  ,  the  estimate  shows  that  under  certain  conditions  on 
the  data,  ajj'  must  grow  quadratically  in  time. 

l'(o,l) 

Theorem  III.  Let  u(x,t)  be  a  solution  of  (l.?4)  with  X(^)  satisfying  (Xl)  - 

(X;J.  Then  if  £(0)  >0  ,  J(u  ,v  )  >0  with  H^u  (x)v  (x)dx  >  J?t{ 0)  (  f5u  (x)dx)1/ 

oo  ^  o  o  'Jo 

-T  oo  —  oo 

we  must  have  on  0<t<tma x 

(2.YJ)  <&V62)  lull'-’  >  4uq)  +23/2ye(o)  v'j(uo)t  +2e(0)t'’ 

l'(o,l) 

i  roof-  We  begin  'with  (L'.'jb),  i.e. 


FT"  -  (v+l)F12  >''oc'(2y+l)F  ,  0<t<t 


max 


where  v  •-  >  0 

iy  cur  hypotneses:  F(0)  >0  ,  hence  3T)>0  s.t. 


F(t,  >0  ,  te  f  0,7)1 
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We  multiply  the  differential  inequality  (2.38)  through  oy 
-Y(F"Y),(t)(F‘^+“)(t))"  ,  t  e  [0,T|] 
and  integrate  over  [0,t]  ,  t<7]  so  as  to  obtain 

(P.51)  f  (F'Y)  '(t)]2  -  Uy2v:?  F‘^’Y+l)(t)  >[  (F"Y)  '(0)]“J  -  4YVf"'  ’Y+l)(0)  >0 


by  virtue  of  the  definition  of  F(t)  and  the  hypothesis  relative  to  the  initial 
data.  Factoring  both  sides  of  (lM  we  have 


(2.52)  [  (F_Y) '  (t)  -2yv  p~(Y+1/fc>;(t)  ][  (F‘Y)  •  (t)  +2Yv  f-(y+iA’)  (t)  3 

>  [(f~Y)’(0)  -2Yv  F'(y+1)(0)][  (f‘Y) '(0)  +?Yv  f-W+1/2)(0)] 


and  thus  as 

(F~y)  '  (t)  =-  Y  F(Y+l)F'(t)  <0  ,  t  e  f  0,rj  , 

(2.53)  (F_Y) '  (t)  < -2Yv  p-(Y+1/2)(t)  >  t  e  10, T]] 

nence,  by  continuity  we  can  not  have  F( T| J  =0  ,  for  any  T1  > 0  .  Thus  F(t)  >0 

0<t<t  and  (2.53)  holds  for  all  t  ,  0<t<t 
—  max  —  —  max 

From  (2.53)  we  obtain  directly  the  estimate 

'2  5M  F(t)  >  (vt  +F1//2(0)  )2  ,  0<t<t 

~  —  max 

and  the  quadratic  growth  estimate  now  follows  from  the  definition  of  F(t) 
and  the  estimate 

l!u]|y  >  ^  [  f  u(y.t)dy)  dx 

.2,  0 

L 


■a 


j2 

3 .  Rlemann  Invariants  and  Finite-Time  Breakdown  of  the  Elec  t_r  ie_Ir.duc  t  ior_  Fit  Id . 

In  this  section  we  otfer  a  brief  demonstration  of  the  fict  that  under  a 

*  1 

slightly  different  set  of  assumptions  on  A  (£),£  £  R  ,  other  than  those  repre¬ 
sented  by  (.XI)  -  (A4),  it  is  possible  in  certain  situations  to  apply  the  Riemann 
Invariant  argument  of  Lax  C 1 3 ]  so  as  to  conclude  that  finite-time  breakdown  of 

u  +—  /ifi'  (u)  u  must  occur,  where  u(x,t)  is  a  solution  of  the  initial- 
t  r  x 

rp 

boundary  value  problem  (1.22);  the  requisite  assumptions  on  >*(£),  however, 

2 

cannot  be  realized  in  the  case  where  A  *  (t )  =  A*  +  A*c  ,  A*  >  0,  A*  ■>  0.  In 
this  latter  situation  it  is  possible,  however,  to  apply  some  recent  results  of 
Klainerman  and  Majda  [21]  as  we  will  indicate  below. 

In  LI 3]  Lax  considers  the  nonlinear1  initial -boundary  value  problem  on 


(3.1) 


ytt(x,t)  =  Khypy^Cx.t), 

y(x,0)  =  y^Cx),  yt<x,0)  =  0;  0  <:  x  s  L 

y(0,t)  -  y(L,t)  =  0,  t  >  0 


This  problem  nay  be  extended  to  a  pure-initial  value  problem  on  R1  *  [0,  °»)  by 
extending  y^(’),  y(’,t)  as  odd  functions  to  (-L,L)  and  then  periodically,  to 
all  of  R  ,  with  period  2L.  By  setting  U  =  y  ,  V  =  y^  the  resulting  extended 
i tit ial -value  problem  on  R*  is  then  easily  seen  to  he  equivalent  to  a  pure 
initial -value  problem  for  a  coupled  quasilinear  system  a.i  RJ  *  i  0,°»)  ,  t.o., 


/u\  /°  -n  r  \ 

V  V  /  ,t  V  -K2  (U)  0  /  V  7  /  .} 


( )  - ; ) 


/  U(x,0)  \  /y0(x)\ 

V  V(x,0)  /  Vo/ 


y^  tlie  extension 
of  yp  to  R^ 


The  eigenvaluer  and  eigenvectors  associated  with  the  system  (3.2) 


are 


9 


respec+ ively , 


n 


( 6fr)) 

•  1  ' 

{  -  =  ±  KCJ)  and 

}  P(b') ) 

.-tk(U). 

irt- 1  thus  the  system  is  hyperbolic  if  and  onlv  if  K"(4)  >  it ,  v  %  <_  R' .  Also,  •»!.« 
system  ;:uy  be  diagonalized  in  a  f  fra  l  iar  way  no  as  to  yield  *h.»  syst •.-.ti 


(3.4) 


V'  +  K(U)'J'  =  n 
V'  -  K(U)U'  =  0 


wl  iere 


3_ 

at 


K(U)  ~  and 
3x 


.  ^  +  K(U) 


D>: 


denote ,  respectively,  ii  1  I'ent. i  it  io- 


along  the  right  and  left-hand  cJ laracterist ios  defined  by  the  onliiury  different  hi  i 
equations  ^  =  ±KU).  dying  (3.4)  one  then  shows  in  the  standard  w«y  that  the 
P.ierrunr.  Invariants 


R(U,V)  =  V  +  /lnJ  K(0  d c 
( .■■ .  5 )  u 

S(U,V)  =  V  -  /Jj  K(?)  dc 

satisfy  R"  =  S'  =  0,  i.e. ,  that  they  aio  coiiStant  along  the  respective  character¬ 
istic  oaves.  It  is  shown  in  [13]  that  with  a  suitable  choice  of  H  =  H(R,A),  th<. 
f  motion  Z  =  exp(A)R^  satisfies  1'  -  -[  (exp(-W) )6,,JZ''  whexe  6R  =  sjp  j'-SW) 
so  that  Z,  and  hence  R  ,  must  breakdown  (blow-up)  in  finite  time  il  3  C  >  *' 
s i .  Ih  that  \  (uxp(-H)6  |  z  C;  this  last  condition,  on  the  other  lurid,  turns  out 
to  he  a  consequence  of  the  assumpt ion  that  lit  >0  :nd.  that 

|DK/an|  >  e  >  d ,  V  U  €  R1 .  rinitc-tim  bre.  dedown  fa-  R  then  Imol  ios  finite-time 
breakdown  for  at  least  one  of  the  second -order  derivatives  y  ,vtt  ct  tie 
solution  y(x,t)  to  the  nonlinear  ini  fial-beun  lory  vain'7'  problem  (,.l),as 


R 

x 


-  R,j  ux  +  Rv 


V 

x 


=  K(1J)U  +  V 
x  x 


-  K(y  )y  t 
x  xx 


xt 
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Suppose  that  we  now  reconsider  the  inriul  mUh  prot  3 .  n,  M 
and  recall  that  as  a  consequence  of  the  fact  that  supp  u  ( x ,  t )  -  (-r, ,rj 


x 

uf  u  (y,t)dy  =  i^(u(x,t)),  =  v1  (u)u  (x,t) , 

_  co  IT  x  * 

•ii(Q)  -  fi*(c)>  t  €  R  .  If  we  set 


for 


x  •  x  where 


(3.6)  v(x,t)  =  /Xco  u^(y,t)dy,  t  i  0 

then,  clearly  ,  v^(x,t)  =  u^(x,t)  and  vt(x,t)  =  /Xqo  ut1_(y,t)dy  =  ip'  (u)uy(x,t) 

Also,  u(x,0)  =  Uq(x),  v(x,0)  =  /Xoo  u+(y,0)dy  =  fXoc  v^tyWy.  Therefore, 
the  initial-boundary  value  problem  (1.22)  for  u(x,t)  is  easily  seen  to  be 
equivalent  to  the  following  initial-boundary  value  problem  for  the  pair 
(u(x,t) ,  v(x,t)) : 


(3.7) 


—  ip'  (u)u 

U  X 


0 


-®  <  X  <  00 

t  a  0 


u(x,0)  =  u^(x),  v(x,0)  -  /Xao  v  (y)dy,  -< *  <  x  <  « 
u(0,t)  =  u(L,t)  :  0,  t  2  0 


The  system  (3.7)  is  clearly  of  the  same  form  as  that  considered  by  Lax  ! 13 J,  i.o 
(3.2),  if  we  assume  that  v  (x)  '•  0,  -«•  <  x  <  «. 


'3 


Thus  in  standard  matrix  form  our  system  { 3.7)  is 

<  x  < 
t  2  O' 


in  comparing  (3.3)  with  (3.2)  we  clearly  have  the  correspondence 
3  1  1 

K  (;)  =  —  if'(0,  S  f  R  ,  and  thus  (3.8)  ->'s  a  hyperbolic  system  if  and  only  if 


ip'CO  =  ?  X  '(0  +  U;)  >0, 


which  is  precisely  hypothesis  (A3).  The  kiemann  Invariants  associated  with  the 
system  (3.8)  are,  clearly,  given  by  the  expressions 


-t(u,v)  =  v  +  — ~  fr,  *V  (p)dp 


(3.9) 


>'  U 


4(u,v)  =  v - —  /q  /PTpTdp 

/  p 


met  they  satisfy  ti'  -  ~  0  along  the  respective  chcu\ietet'i  sties  given  1  y 

/^r  :•)  .  .  _  j  fy  ) 


-  ±  /PTuT  where 
v/  ~  Ij- 


St  /  1 


I  cX 


and 


St 


•  •  /  the  result t  in  [13], which  we  have  described  dhow .  f  su.  .e-t  i  v;  breakdown 
(blow-up)  of 

K  -  V  +  - -  v  lj/'(li)u 

x  x  r  x 

/U 


-•  Ut  + 


—  y~ i 1 (u  >u 
4i 


(3.13) 
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w  if,  V  ;  f  R  (this  is  actually  ate  cit'd  only  7  ^  sufficiently, 

l(  L  Ik-  Initial  data  is  small,  as  we  indicate  below), 


v"(Q 


=  — 1 - •  wet)  !  >  c“ 

?Suip'(ZJ 

for  some  e  >  C.  Using  the  relationship  between  X(')  and  \)  this  la:  t 
condition  is  equivalent  to  the  requirement  that  satisfy,  for 

e  >  U, 


(AS) 


+  2X’(r,)|  >  e/;X'(c)  +  A (O, 


It  also  follows  from  the  work  of  Lax  [13]  that 


t  *  4 

max  max  Uq  ( x )  ip"  CO) 

A  A  •% 

As  (Jt'(c)  =  tt'(0  +  A(C>  aTd  -  cA"'0  +  2 A’(OsVC  «  R  ,  we  clearly 

A 

must  require  that  AC?)  satisfy 


(Ab) 


0  <  A ( 0)  <  «-,  C  <  A '  (U>  <  «,  j  A"(0)  I,  < 


in  which  case  t  ^ - .2j5L„. 

tiax  max  u;,(x> 


Now ,  i 


\J^l  \ 

0W“ 

tnat  the  a  priori  estimate  b(x,«  (x.  t  j  ..  *..t> 


x 

0:1  (-«<,*)  x  fO,tmax],  By  virtue  of  the  definitions 

e-itimates  held  for  u(x,tj,  v(x,t). 


a  simple  ratter  to 

> '  x , 0 )  ,  -*■  sup  li  <x,n)  | 
x 

1.9)  ’Similar  a  priori 


J7 


and  thus  j  u(x,  t)  1  and  '  v  <  y ,  t )  ,  will  be  sma  1 1 

if  sup!u(x,0)|  is  sufficiently  small.  It  then  follows  that  (/  5)  need 

only  hold  for  some  e:  0  ana  all  \  (  R  *  which  are  ;uf  f  ir  i«.*nt  ly  small 

*  *  *  *  * 

in  magnitude.  Now,  if  >.  (s)  =  aq  +  ']'-•>  aq  >  11 »  ■'  ’  *■',  '  '  9*  then 

(X5)  becomes 

*  /*  * 

(3.11)  2>,  -  i  /2‘t,  r\Q 

which  is  certainly  satisfied  for  | C !  A.  A  >  0,  if  t  ■  0  is  chosen 

*  *  2 

sufficiently  small.  In  this  case  ’1- ( r )  =  +  a  r  so  that  v  (0)  = 

k 

2X^  ^  0  and  the  system  (3.8)  is  genuinely  nonlinear.  However,  if 

*  -k  k  2  *  * 

X(-r.)  =  Aq  +  X  t,  ,  A  >0,  a  0  then  (a 5)  becomes 

(3.12)  6*2k|  2  e/3X2k  +  >*0 


which  cannot  be  satisfied  for  any  f  >  0  even  as  U(  0.  The  essential 

problem  here  is  the  lack  of  genuine  nonlinearity  vis  a  vis  ( ' )  which  in 

*  *  3  A 

this  case  is  ip(rj  =  1  ;  +  >  ^  (so  that  i'  (O  =  6‘0C  which  vanishes  at 
i,  =  0).  In  order  to  obtain  finite-time  breakdown  of  smooth  solutions  to 
initial-value  problems  for  the  quasiline', r  system  (  i.7)  on  ,  when  tin1 


data 

u(x,0) 

,  v(x 

,0) 

has  ccmnact  support  in 

r\ 

and  i' ( )  -  sl(  \)  = 

V  + 

?  r3 

2  ’ 

* 

A  0  > 

o, 

>  0,  we  may  appeal 

to  3 

t  ‘cent  result  of  S. 

Klainerman  and  A.  Majda  [21;.  Stater,  in  terms  of  the  in •  t > a  1-va lue  problem 
for  the  diagonal  system  associated  wifi:  me  Kicmann  ianr.3,  i.e.. 


(  l.  13) 


r  *■ 


d 1  A’  A T 
dt  /  \i  Ax 


d±  /£_  U 

at  /  p  7 


4' 


0 
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where  K,i>  are  given  by  (3.9),  the  relevant  result  in  f  21  ]  says  that  any 

Cl  solution  of  (3.13)  with  C1  inital  data  4  (x)  =  T(x,0),  6  Q  ( x )  =  4(x,0), 

having  compact  support  in  r\  must  develop  singularities  in  the  first 

derivatives  A.  ,  6^  in  finite  time  provided  (•'((;)  is  not  constant  on  any 

open  interval;  this  last  condition  is,  of  course,  equivalent  lo  having 

'('"(C)  ^  0  on  every  open  interval  -  which  is  certainly  true  for  the  cubic 

<KO  of  our  example.  It  is  worth  noting  that  as  with  the  work  of  Lax  [li  l, 

the  work  of  Klainerman  and  Majda  [21]  was  done  with  the  nonlinear  wave  equation 

in  (3.1)  in  mind  and,  as  in  [13],  the  development  of  singularities  in  the 

first  derivatives  A  ,  4  leads  to  tne  prediction  that  solutions  of  the 

x  x 

initial-value  problem  for  (3.1)  with  C  initial  data  having  compact  support 

in  R1  must  develop  singularities  in  the  second  derivatives  y  ,  y 

xx  xt 

However,  as  is  the  case  with  the  work  in  [13],  the  results  in  [21]  now  pre¬ 
dict  that  solutions  of  the  initial-value  problem  (1.22),  with  C1  initial-data 
having  compact  support  in  R1,  must  develop  singularities  in  the  first  deriva¬ 
tives  u  ,  u  . 


Remarks .  In  closing  we  offer  a  few  comments  concerning  the  problem  of  proving 
global  nonexistence  of  smooth  electric  induction  fields  of  the  form  (1.20)  in 
a  finite  rod  occupying  the  configuration  0  x^  L.  The  relevant  one¬ 
dimensional  equation  is  still  (i.2L)  but  now  we  must  take  account  of  the  impli¬ 
cations  of  the  boundary  conditions  (1.16),  assuming  as  before,  that  the  rod  is 
embedded  in  a  perfect  conductor.  At  the  planar  boundary  at  x,  =  0, 
v  =  (-1,0,0),  £  ■*  (0,1,0)  and  thus  by  (1.16a)  with  =  0,  f(x2*x^)  - 


(3.14)  D(x,t)  •  v 


xr° 

t>0 


(0,D  (*  ,t),0)  •  (-1,0, 0)1 


V° 

t>0 


.so  t'lat  (x,t),  the  surface  charge  density  at  y  =  0  rwi-:i  vanish  t  ;.r  : 

t  >  0.  An  analogous  result  holds  at  x,  -  L  where  y  =  ( j ,0,0; .  In  order 
to  satisfy  the  boundary  condition  (1.16d)  along  the  planar  face  at  .  =  0, 
for  t  >  0,  we  require  that 

(3.15)  X(D(x,t))D(x,t)  •  tl  n 

. . h=0 

1  t  0 

*  (0.?(D2(*1.t)D2(x1.t).0)  •  (0,l,0)jx=0 

■t>0 

=  *(D2(x1,t))D2(x1,t)|  0 


from  which  it  follows  that  D2(0,t;  =0,  t  >  0.  In  an  analogous  manner  we 
have  D2(L,t)  =  0.  In  place  of  the  initial-value  problem  (1.22)  for 
u(x,t)  D^Cx^.t)  we  then  have  the  initial-boundary  value  problem 

2  2 

u  ~  CuX(ti)),  (x.t)  '  i  0 . :  ••  x  10, T). 

St '  Sx*” 

<  , 

u*.0,t)  =  u(L,t)  =  0,  t  -0 

n(x,0)  =  u^fx),  (x,0)  =  vq  1 K 1  ’  c  x  •  !. 


and,  in  addition,  because  of  the  embedding  of  c.oe  rod  in  a  perfect  conductor, 
supp  u  c  r0,L].  The  principal  difficulty  that  arises  in  trying  to  apply  the 
analysis  of  both  this  and  the  previous  section  to  ei'iar  (3. Id)  or  the  m<n‘ 
initial  value  problem  that  results  by  making  the  asn.1  e:.tcr.si-,  as  of  the  initi 
data,  first  to  [-L.L]  and  then  to  ail  of  R*  with,  period  ?L,  revolves 
around  dealing  with  the  integral  u  fv,t)’.y  which,  in  the  analysis  ot  th 


infinite  rod,  is  equal  to  —  i)/'(u)u  (x,t).  If,  as  is  customary  in  trying  ;<■ 

U  X 

prove  breakdown  of  smooth  solutions,  we  assume  that  u(-,t!  is  of  class 
2  1 

C  ,  for  all  t  >  0,  on  R  then  integration  across  the  planar  undnry  at 
x^  =  0  forces  upon  the  analysis  the  a  priori  assumption  that  not  on:y 
u(0,t)  =0,  t  >  0,  but  also  u  (0,t)  =  u  (0,t)  =0,  t  >  0;  is  un- 

X  XX 

likely  that  any  classical  solution  of  (3.16)  could  exist  under  sure  ci^cum- 
(4) 

stances.  We  hope  to  address  the  problem  of  global  nonexistence  oi  smooth 
solutions  to  the  initial-boundary  value  problem  (3.16)  in  a  future  paper. 


w 


The  author  is  indebted  to  Prof.  Morton  Gurtin  for  this  observation. 
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